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Abstract 

A completely antisymmetrized Green's function approach to the inclusive quasielas- 
tic (e, e') scattering, including a realistic one-body density, is presented. The single 
particle Green's function is expanded in terms of the eigenfunctions of the nonher- 
mitian optical potential. This allows one to treat final state interactions consistently 
in the inclusive and in the exclusive reactions. Nuclear correlations are included in 
the one-body density. Numerical results for the response functions of ^^O and ^'^Ca 
are presented and discussed. 
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1 Introduction 



The one-body mechanism gives a natural interpretation of the inclusive elec- 
tron scattering in the quasielastic region. However, in order to explain the ex- 
perimental data of the separated longitudinal and transverse responses more 
complicated mechanisms are needed. A review of the experimental data and 
their possible explanations can be found in Ref. [1]. Thereafter, only a few 
experimental papers were published [2,3], while new experiments with high 
resolution are planned at JLab [4]. 

Many papers were published in order to explain the problems raised by the 
separation, i.e., the apparent lack of strength in the longitudinal response and 
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the apparent excess of strength in the transverse one. Among them, the more 
recent ones are concerned with the contribution to the inclusive cross section 
of meson exchange currents and isobar excitations [5-8], with the effect of 
correlations [9,10], and the use of a relativistic framework in the calculations 
[7,11]. 

At present, however, a consistent and simultaneous description of the longi- 
tudinal and transverse response functions is not available. A possible solution 
could be the combined effect of two-body currents and tensor correlations 
[12,9,13]. 

A peculiar problem of inclusive electron scattering is the treatment of final 
state interactions (FSI), since they are essential to explain the exclusive re- 
action with one-nucleon emission, which is the dominant process contributing 
to the inclusive reaction in the quasielastic region. The large absorption due, 
e.g., to the imaginary part of the optical potential, produces a loss of flux that 
is appropriate for the exclusive process, but inconsistent for the inclusive one, 
where the total flux must be conserved. This conservation is preserved in the 
Green's function approach considered here. This result was originally derived 
by arguments based on the multiple scattering theory [14] and then by means 
of the Feshbach projection operator formahsm [15-18,11]. 

The spectral representation of the single-particle (s.p.) Green's function, based 
on a biorthogonal expansion in terms of the eigenfunctions of the nonhermitian 
optical potential, allows one to perform explicit calculations and to treat FSI 
consistently in the inchisive and in the exclusive reactions. In previous papers 
[17,11] the Green's function approach was used both in a nonrelativistic and 
in a relativistic framework to perform explicit calculations of the inclusive 
response functions. 

Two issues are the main goal of this paper: a completely antisymmetrized 
presentation of the Green's function approach, and the effect of nuclear corre- 
lations, which are included in the model through the use of a realistic one-body 
density matrix (ODM). In the previous application of the method [17,11] cor- 
relations were neglected. The response functions were given by a sum over 
the residual nucleus states restricted to be s.p. one-hole states in the target. 
A pure shell model was assumed for the nuclear structure and therefore the 
overlap functions between the target and the residual nucleus were given by 
phenomenological s.p. wave functions with a shell-model spectroscopic factor. 
In the present work we are able to include partial occupation numbers, and 
therefore correlations, through the natural expansion of the ODM. 

The deflnitions and the main properties of the quantities involved in the model 
are given in Section 2. In Sec. 3 the antisymmetrized Green's function approach 
is developed and the inclusive cross section is expressed in terms of the ODM. 
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In Sec. 4 the Green's function is calculated in terms of the spectral represen- 
tation related to the optical potential. In Sec. 5 some models describing the 
ODM including correlations are briefly reviewed. In Sec. 6 the results of calcu- 
lations for ^^O and ^'^Ca are reported and compared with some experimental 
data. Summary and conclusions are drawn in Sec. 7. 



2 Definitions and main properties 

In this Section we collect most of the definitions which will be useful in the 
treatment of the inclusive (e, e') reaction, using the Green's functions approach 
along the lines of Refs. [17,18,11]. 

2.1 One-body density, Green's functions and related quantities 

We deal in this paper with three different systems: the Z-proton residual 
nucleus, the {Z + l)-proton target nucleus, and the [Z — l)-proton system 
obtained removing a proton from the residual nucleus. The contribution of 
the neutrons, here disregarded in order to simplify the formalism, can be 
introduced in an obvious way. Let T and H be the kinetic energy and the 
Hamiltonian operators acting in the Hilbert spaces 7Y^, TY^"*"^, and Ti.^~^, cor- 
responding to the three given systems and written in the Fock's formalism. 
Here we are interested only in the following eigenvectors of H: i) the eigenvec- 
tors |n), representing the bound states of the residual nucleus with energy e„, 
ii) the eigenvector \ipo), representing the initial state, which is the ground state 
of the target nucleus with energy Eq, and iii) the eigenvectors lip/), represent- 
ing all the final (bound or unbound) states of the target nucleus, including 
l^/'o). These eigenvectors are properly orthonormalized. Relatedly, we define 
the hole overlaps, referred to \ipo), 

{r\(f)n) = (nlorlV'o), (1) 

and the particle overlaps, referred to |n), 

{r\Xf,n) = {nWrli/jf), (2) 

where annihilates a proton at the point r. The corresponding annihilation 
operator of a proton of momentum p is 

ap = (27r)"^/^ j dr exp(-ip ■ r) ar. (3) 

For sake of simplicity, we consider spinless protons, we assume that \ipo) and 
|n) are not degenerate, and we omit any degeneracy index in the unbound 
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states \ipf)- When referring to quantities related to n — 0, the index will be 
usually suppressed in the following Sections. 

2.1.1 The one-body density matrix 

The one-body density matrix related to \ipo) is defined as 

{r\K\r') = (V'ol^arlV'o). (4) 

It is real and symmetric by the exchange r < — > r', since lipo) is not degenerate. 
The corresponding operator K in L^(R^) is nonnegative, self-adjoint, and has 
a finite trace equal to Z+l. The eigenvalue equation 

K\u^) = ny\ui,) (5) 

defines the natural orbitals Uu{r) — {t\ui,) and the related occupation numbers 
nj/, which satisfy the relations 

< 1, = Z + (6) 

V 

The operator K is invertible in the subspace of L^(R"^) spanned by the orbitals 
Uy^ with ^ 0. If = is not an eigenvalue of it is necessarily an 
accumulation point of eigenvalues. Thus K is fully invertible, but is a 
rather complicate operator. 

2.1.2, The one-body Greenes junction 

At complex energies ^, the one-body Green's function related to \n) is the sum 
of the particle and hole Green's functions: 

(r|^„(z)|r') = {t\Q^\z)\t') + (r|^W(z)|r') (7) 

with 

kT\G^^\z)\T') = {n\ariz-H + e^)-'al,\n), (8) 

(rl^W(z)lr') = {n\al,{z + H- er.)-'ar\n). (9) 

These Green's functions are symmetric by the exchange r < — > r', since \n) is 
not degenerate. 

At real energies E, we consider the retarded Green's function instead of the 
time-ordered one, since it is more convenient in the formalism developed below. 
It is defined as 

QniE) = lim GniE + i'q). (10) 

In this equation, as well as in all the next equations involving rj, the limit is 
understood in the weak sense, i.e., it must be performed after inclusion into a 
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scalar product between normalizable vectors. Henceforth, the symbol of limit 
will be usually understood. 



2.1.3 The spectral function 

The Green's function Qni^) is analytic in the complex plane except for cuts 
and poles on the real axis. Hence, it can be written as 



where 



'^"^^^ = ^^[Gn{E-iv)-Gn{E + iv)] = ^^[gi{E) - Qr^iE)] (12) 
is the spectral function, which satisfies the sum rule 

/+0O 
dESn{E) = 1. (13) 
-oo 

Since Sn{E) is nonnegative, Eq. (11) and (13) show that 

g^izM^o ^ |(/.) = o (14) 

for z complex. Hence, Qniz) is fully invertible. 
2.1.4 The self- energy 

At complex energies z, the Green's function can be written as 

where T is the one-body kinetic energy and A1„ the one-body self-energy. The 
related Hamiltonian hniz) is defined as 



hn{z) =T + Mn{z) ^Z- [Sn{z)]-\ (16) 

Note that Qn{z) is fully invertible, so that M.n{z) has no mathematical draw- 
backs related to undue restrictions of its domain. In contrast, QI^\z) is only 
partially invertible if 1 is an eigenvalue of the density matrix associated with 
\n). This produces a restriction of the domain of the related self-energy M.^\z) 
which may lead, e.g., to an incorrect Dyson equation. The same drawback af- 
fects M'il'Xz) if is an eigenvalue of the density matrix. 

The analyticity of Qn{z) induces similar analyticity properties into hn{z) and 
M.n{z). They are analytic in the complex plane except for cuts and poles on 
the real axis, which are different from those of Qn{z). 
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The one-body self-energy at real energies M.n{E) is defined as 

Mn{E) = Urn Mn{E + i7]) (17) 
»?— >+0 

In the Appendix B of [21] it is proved that Gn{E) is fully invertible, except 
for the values of E corresponding to the poles of A4n{E), and that holds the 
relation 

The related Hamiltonian hn{E) is defined as 

K{E) = T + Mn{E) = E- [gn{E)]-\ (19) 



2.2 Self-energy and extended projection operators 



In his pioneering papers [19,20] Feshbach introduced a projection operator 
formalism to obtain a closed expression for the theoretical optical-model po- 
tential. Here, we use an analogous procedure for the self-energy, based on 
different projection operators. More details can be found in Ref. [21]. 

Let us introduce the vectors aj.\n), where 

ar = a-r + al . (20) 

They belong to the Hilbert space Ti^'^"*"^-* © 7i^^~^^ and form an orthonormal 
set. In fact one has 

{n\arar'\n) — {n\aral, + al,ar\n) — S{r — r'), (21) 

where we have used the property that (n|ajar'|n) is symmetric since \n) is not 
degenerate. Therefore the operators 

Pn^ j dr ar\n) {n\ar ^ j dpap\n){n\ap , Qn^l-Pn (22) 

are projection operators in 1-6'^^^'' © T-C^^~^\ Such operators, which will have 
an important role in Sec. 3, are called "extended projection operators" in order 
to distinguish them from the Feshbach's ones. For sake of simphcity, we define 
here P„ in the simpler case of spinless particles. The necessary changes to 
include spin and isospin variables can be found in the Appendix A of Ref. [21]. 

One can easily check that the correspondence 

Pr,\(f>) ^ {n\ar\cf>), (23) 

PnOPn < — ^ {n\arOar'\n), (24) 
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where |0) and O arc vectors and operators in Ti.^^~^^^ © defines an 

isomorphism between the Hilbert space of the vectors Pn|0) and L^(R"^). This 
means that every property involving vectors and operators of one space implies 
that the same property holds in the other space. 

The above defined isomorphism is useful to develop a more synthetic formal- 
ism for the quantities defined in Sec. 2.1. To this extent, we introduce the 
Hamiltonian-type operator defined by the relations 

Hn = H-en in Ti^^+i) , = en - H in n^^-'\ (25) 

and extended to the whole space 7i^^~^^^ © H^^~^^ by linearity. Moreover, we 
introduce the related quantities 



Z - -tin 

Sn{E)^5{E - Hn) = ^{Gi{E) - Gn{E)). (26) 
By the definition oi S{E — Hn) one has 

/+00 
dESniE)^!. (27) 
-oo 



Due to the orthogonality between vectors related to different numbers of pro- 
tons and due to the symmetry properties of the particle and hole components 
of Gn{z), Eq. (7) is reduced to the simple expression 

{r\gn{z)\r') = {n\arGn{z)ar'\n). (28) 

Thus, Qn{z) is the one-body operator isomorphically corresponding to the 
many-body operator PnGn{z)Pn- Due to Eqs. (12) and (26) it follows that 
Sn{E) corresponds to PnS{E — Hn)Pn- Therefore, one has 

{r\Sn{E)\r') = {n\a,5{E - Hn)a,'\n). (29) 

As the handling of equations is often simpler when one deals with many-body 
operators of this type, we introduce, in keeping with the definitions of Sec. 2.1, 
the many-body definitions of the Green's function, the spectral function and 
the self-energy: 
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Gniz)^PnGn{z)Pr, , Gn{E)^P^Gr,{E)P^, (30) 

S^{E) = PJ4E)P^ = P^S{E - H^)P^ = :^[Gi{E) - G„(E)], (31) 

K{Z) =Pn{z- [Gn{z)]-') Pn , K{E) = P„ [E - [G'n(^)]-^) Pn , (32) 

M^{z)^hn{z)-Pnf^Pn , M^{E)^hn{E)-P^f^P^, (33) 
where 

f 

Tn = J dfeakin) — (n|ak, (34) 

is the many-body operator which corresponds isomorphically to the one-body 
kinetic energy T. 

The closed expression of the many-body self-energy is deduced in Sees. 3.1-3.3 
of Ref. [21] and is: 

Mn{z) = Pn[Hn ' fn + HnQn{z - Q nH nQ nY^ Q nH n]Pn • (35) 

For a two-body potential 

V = ^ J drds dr' ds' alalV{r, s, r', s^Og'ttr', (36) 

the corresponding one-body self-energy Mniz) is the sum of the static and 
dynamical parts 

(r|M(fV) = (^l^raj:, + alJr\n), (37) 

(r|Mf V) = {n\{Jr + Jl)Qn{z - QnHnQnY^QniJr' + 4)1^), (38) 

with 

Jr = [ar,1/]. (39) 

The structure of M.n{z) is identical to that of the Feshbach's potential, but 
its one-body expression is different. For instance, the Feshbach's potential is 
not symmetric for the exchange r < — > r' . In contrast, A4.n{z) is symmetric, 
as one can check working out the right-hand side of Eqs. (37) and (38). 



3 Green's function approach to inclusive (e, e') reactions 

3.1 Hadronic tensor 

In order to avoid complications of minor interest in the present context, we 
omit recoil corrections, use the one-photon exchange approximation with one- 
body currents, and consider only spinless point-like protons. 
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The inclusive cross section for the quasielastic (e, e') scattering on a nucleus 
of Z + 1 protons is given by 

a = K(2eLi?L + i?T), (40) 
where is a kinematical factor and 

measures the polarization of the virtual photon. In Eq. (41), 9 is the scattering 
angle of the electron and = c<j^ — where {uj.q) is the four-momentum 
transfer. All nuclear structure information is contained in the longitudinal and 
transverse response functions, defined by 

R^{uj,q) = W^\uj,q) + W^\uj,q) (42) 
in terms of the diagonal components of the hadronic tensor 



W^'^iu, q) = {^^\r\q)5{u + Eo- H)r{q)\i;,) 

= ±{MJ'\<l)\^s){i^}\nQ)\^o)8{u + Eo- Ef). (43) 

Here J^{q) is the one-body nuclear charge-current operator, I'i/'o) is the ini- 
tial state of the {Z + l)-protons target nucleus of energy Eq and |^/) is the 
corresponding final state of energy Ef = Eq + uj. The target ground state 
\iPq) is assumed to be nondegenerate. The sum runs over the target bound 
states and the scattering states corresponding to a proton scattered from the 
residual nucleus in a bound state \n) or in an unbound state \e). For sake of 
simplicity, the degeneracy indices are suppressed. All the states are properly 
antisymmetrized. Accordingly, the nuclear Hamiltonian H and the current 
components J^{q) are understood as second quantization operators. 

3.2 Energy sum rules for the incoherent and the coherent contributions 

Acting on l-^o) the scalar component 

J°(g) = J dr exp(iq ■ r) ajor (44) 

yields the wave function 

z+i 

Jii^^ ri)ipo{ri, rz+i) , J-iq, Vi) = exp(ig • r^). (45) 

i=l 
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Similar expressions are obtained for the other components J'^ of the CTirrcnt 
operator. The hadronic tensor of Eq. (43) is the sum of the incoherent contri- 
bution 

W^Ziu, q)^iZ + l){,l;o\J^\q)5iu; + E, - H)J^{q)\i^,) (46) 
and the coherent one 

2+1 

W^cTh(^, q) = (Z + 1) ^ {MJ'i\q)K^ + Eo- H)Jt{qMo). (47) 

i=2 

The non-energy-weighted electromagnetic sum rule deals with the integrated 
strength 

W^(q) = r du;W^^iu;,q) = {,l;o\J'\q)nqMo) , (48) 
J-o 

where the lower integration limit means that the integral includes the S- 
singularity at cu—O (full sum rule inclusive of the elastic contribution). In 
the scalar case one has 

^^M = Z + 1, (49) 

z+1 

^ooM = iz + i) T.iHJi\Q)J'i<im. (50) 

i=2 

In the Fermi gas model Sp°ij(q) is negligible for |q| > 2A;f, where kp is the 
Fermi momentum. This is still approximately true in presence of correlations 
[Refs. [22,23]]. Therefore one has 

S°°(q)~(Z+l) for \q\>2kF 

3.3 Projection operator method 

The approach developed in this Section differs form previous treatments [15,17,18] 
for three reasons. 

i) We want to include the effects of the final state interaction in terms of the 
self-energy, rather than in terms of the Feshbach optical potential or of the 
particle self-energy. In fact, the full self-energy is a more fundamental quantity, 
has no mathematical drawbacks, and is more closely related to the empirical 
optical-model potentials. Therefore, we shall not use the Feshbach's projection 
operators as in Ref. [17], but the extended ones of Eq. (22), already used in 
Ref. [18]. 

ii) We want to include the interference between different channels \n), ne- 
glected in Ref. [18], which is useful to express the hadronic tensor in terms of 



(51) 
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the local potential equivalent to the self-energy. This requires some changes 
in the approach of Rcf. [18]. which give in practice the same result when the 
interference effects are disregarded. 

iii) We want to introduce different approximations for the elastic and inelastic 
contributions to the hadronic tensor, when dealing with the dependence on 
the state of the residual nucleus. 

In order to keep the treatment as simple as possible, we shall only refer to the 
scalar component VF°°(a;,qf), with J°(q) written in the momentum represen- 
tation, i.e. 

J°(q) = J dpatap_q. (52) 

Disregarding the upper indices of the hadronic tensor and inserting the com- 
pleteness relation of the residual-nucleus states, Eq. (43) yields 

W{u;,q)^i2 [ dp{,Po\al-^\n){n\apS{u; + Eo-H)j\q)\'il;o), (53) 

where the sum is understood over \n) and |e). Equation (53) was obtained in 
Ref. [18] in a more complicated way. Here, the same result is recovered with the 
insertion of a completeness relation into the second quantization expression 
of J^{q) and has the same generality, i.e. can be applied to every one-body 
charge- current operator. 

As was done in previous papers, we shall now disregard the contribution of 

the continuum states of the residual nucleus. Although this approximation is 
correct at the energy and momentum transfers considered here, we shall re- 
cover the continuum contribution in Sec. 3.4.3. Equation (53) is approximated 
with 

W{u,q)^Y.^e f dp(V'o|4-q|n)Hap% + £;o-//)J°(q)|V'o), (54) 

where the real part has been extracted in order to restore the real nature 
of Wiuo^q)^ which can be lost after the approximation. This is equivalent to 
introduce the truncated completeness symmetrically in J^\q) and in J^{q). 

Now, the treatment of Ref. [18] is slightly modified in order to obtain expres- 
sions directly related to the self-energy, according to the above items i) and 
ii). This requires the introduction of the vectors (a^ + ap)|n), instead of a,p|n), 

and the use of Eq. (25) in order to feature Hn- To this extent, we add to the 
second factor of Eq. (54), written as 

(n|ap 5{u + Eo-en- ^„) J°(q)|^o), (55) 
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the term 

{n\al8{uj + Eo-tn-Hn)J\q)\i'o) = {n\al5{uj + Eo-2en+H)j\q)\ij^), (56) 

which is nuU, since it is a scalar product between states with a different number 
of protons. Moreover, we add to the terms (55) and (56) the corresponding 
ones with ui replaced by —a;, i.e. 

{n\a^5{-uj + Eo-en-Hn)J\q)\il^^) = {n\a^5{-uj + E^- H)j\q)\'il^o). (57) 

which is null for a; > because the 5-function is not fed with the eigenfunctions 
of H, and 

(n|aj, 5{-uj + Eo-€n- ^„) J°(qf)|V'o) = 

{n\al5(-u; + Eo - 2e„ + //) J°(q)|Vo), (58) 

which is null for the same reason as the term (56). The terms (57) and (58) 
are necessary to obtain a null sum rule for the interference contribution, as 
will be apparent below. We emphasize that the term (57) can be added only 
for ujy^O. 

After addition of the null terms (56)-(58), Eq. (54) reads 

W{u;,q)=A{u;,q)+A{-u;,q) (59) 

with 



"p-q 



A{uj,q)^Y.^^ dp(V'o|4- 

X (n|(ap + at.) S{uj + - - Hn)Aq)\A)- (60) 

In Eq. (59), as weh as in the following analogous expressions, it is understood 
that the ^-singularity at = must be retained only in A{u!, q) to avoid a 
double counting. Therefore, one has 

/oo /■oo 
dujW{uj,q)= dujA{uj,q). (61) 
-0 ^-oo 

Now we use the extended projection operators, defined in Eq. (22), 

Pn = j dpQ;p|n)(n|Q;p and = 1 - P„ (62) 

to separate in Eq. (60) a direct {AP) and an interference (^4^) term, such that 

A{uj,q) = A^{uj,q) + A\uj,q) (63) 
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with 



A''{u;,q) = J2Re [ dp 



n) 



X {n\{apPn5{u; + - - Hn)PnJ\q)\i^o), 



(64) 



^^(^,9) = ZlR'e / dp(V'o|ap-q|ri) 

n 



X 



(n|(apP„5(a; + - - Hr,)QnJ\q)\^o). (65) 



Accordingly, the hadronic tensor W{ijj,q) is decomposed into a direct and an 
interference contribution: 



W'iuj, q) = A^iu, q) + A^{-u, q) (66) 

W\u,q)^ A\uj,q)+ A\-u,q) (67) 

We remark that does not contribute to the sum rule, due to the relation 
PnQn — ^ and to the relation 

/oo ^ f-OO 

duj5{uj + Eo-en-H^)= dE 6{E - H^) = 1, (68) 
-oo J —oo 

which yields 

/CO POO 
du}W\u;,q)= duo A\u;, q) = 0. (69) 
-0 ^-oo 

Note that the terms (57) and (58), which yield the contribution of the negative 
energies to the second integral, are essential to fulfill Eq. (69). 

We emphasize that the insertion of the terms (56)-(58) does not change the ex- 
pression of W{uJ, q), but produces an effect on its decomposition into W^{uj, q) 
and W^{u!,q). As a consequence, the present decomposition is slightly differ- 
ent form that of Ref. [18]. The term (56), absent in Ref. [18], modifies the 
decomposition at very low energies, i.e. at a; < 2e„ — Eq — eg"^, where Cg"^ 
is the ground state energy of the {Z — l)-body Hamiltonian. The term (57) 
produces no effects. The term (58), which is also present in Ref. [18], differs 
here by a shift of 2{Eo — e„) in the value of the argument. This shift is the only 
difference from the treatment of Ref. [18] at the energies u; considered in this 
paper. The changes introduced here influence only the separation between the 
direct and the interference contributions, allowing a more appropriate treat- 
ment of the latter, which in Ref. [18] is totally neglected. The effects of these 
changes on W^{LJ,q) are in practice negligible at high momentum transfers, 
as it will be seen in Sec. 3.4.5. 
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3-4 Direct contribution to the hadronic tensor 



3.4-1 One-body expression of the incoherent and coherent contributions 

Using (62), Eq. (64) is expressed in terms of vectors and operators of a s.p. 
Hilbert space as 

A^{u, q)^Yl M<l>n\fHQ)SnM\^n{q)) , = + £^0 - en, (70) 

n 

where |0„) and iS„ represent the hole overlaps and the one-body spectral 
functions, defined by Eqs. (1) and (29), respectively, and 

{P\f{q)\<pn) = (P - q|0n) = (n|ap_q|V'o), (71) 

{p\^niq)) = {nKJ'iqMo). (72) 

Using the relation 

[ap,J°(q)]=ap_q (73) 

in Eq. (70), the contribution to |$„(q)) due to the many-body current J^{q) is 
split into the contribution of a single proton and that of all the other residual 
Z protons, as 

{p\Mq)) = iPlfiQMn) + {p\S^n{q)), (74) 

where 

(p|(5$„(q)) = {n\.f{q)a^\i;o) 
. z+i 
= VZ+1 / dp2...dpz+i{n\p2...Pz+i){p,P2-Pz+i\ J2 Jii^Mo)- (75) 

•' 1=2 

From Eq. (74), we express W^{u!, q) as the sum of the terms 
Wi^,(u;, q) = Af,,{u;, q) + Al,{-u, q), 

^lA^. 9) = E M<t>n\f{q)SnMf{q)\<l>n) (76) 

n 

and 

W^Dja;, q) = A^Jo;, q) + A^,, (-c., g), 

>l?oh(^,q) = ERe(0n|/^(g)5n(^n)|<^*n(g)). (77) 

n 

The decomposition (74) is exactly equivalent to that of Eqs. (46) and (47). In 
keeping with this correspondence, we call and H^coh "incoherent" and 
"coherent" contributions to . 
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3.4-2 Dependence on the state of the residual nucleus 

No information is available concerning the spectral functions Sn{E) related to 
the excited states of the residual nucleus. Therefore, we are obliged to express 
Sn{E) in terms of So{E). According to previous papers [24-26], we assume: 
i) Sn{E) only differs from So{E) by an energy shift 5e„, ii) the contributions 
of the various shifts can be taken into account by means of a single average 
shift e. Thus, in Eqs. (70), (76), and (77) we set 

Sn(u; + Eo- en) ^ So(u) , u ^ u + Eo - eo - e. (78) 

The choice of e will be discussed later. 



3.4-3 Inclusion of the unbound states of the residual nucleus 

For a single state |e) of the continuous spectrum, one can define neither an 
extended projection operator nor a Feshbach's one. This is the reason why 
the states |e) are usually neglected in the treatments based on the Green's 
function approach. Even if the continuum contribution is negligible at the 
energy and momentum transfers considered here, the lack of completeness 
resulting if one considers only the bound states \n) appears unsatisfactory 
from a conceptual point of view, mainly because the sum rule is not exactly 
fulfilled. This drawback can be easily eliminated if Eq. (78) is used. In fact, 
from Eq. (70), by expliciting the scalar products and using Eqs. (71) and (72), 
we have 

A°(a;,qf) = $^Re / dpdp'(V'o|aJ,_q|n)(n|ap, J°(q)|V'o)(p|5o(a;)|p'), 

uj — u} + EQ — tQ — e- (79) 



Here, it is quite natural to recover the contribution of the continuous spectrum 
after making the substitution 

\n){n\ — > + / de |e)(e| = 1 (80) 

n n 

which yields 



^°(cu,g) = Re j dpdp'(7Ao|ai-qV'^°(9)lV'o)(p|«So(cD)|p') 

^ReTx[K^S^{Co)l (81) 



where 
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(p'l ifqlp) = ('0o|ap-qV>^°(9)lV'o) = J dfc('0o|aJ,_qap/4ak-q|V'o>- (82) 
Analogously, using Eq. (73), one has 

Af^^{u,q)^ReTr[K^So{Co)] (83) 

with 

{P'\K^\P) = (^0|4-qV-ql^0) (84) 

and 

A°,(u;,q) = ReTr[AK^Som (85) 

with 

(p'lAXqlp) = j dfc(V'o|ap-q4«k-qVlV'o)- (86) 

Since -fTq and Sq are nonnegative (hermitian) operators and i^q is self-adjoint, 
Tr[A'qiSo(ti))] is nonnegative, as it follows easily by expressing the trace on the 
basis of the eigenvectors of K^. Thus one has 

AZ^{u;,q)^T:i[K^So{u;)]>0. (87) 

The expression of 

W'^iu, q) = ^°(a;, q) + A°(-a;, q), (88) 

obtained from Eq. (81), and the corresponding expressions of W^^^{uj,q) and 
W^Y^{uj, q), are essentially the same as those obtained in Sees. 9.1.1 and 9.1.3 
of Ref. [18]. This can be easily checked noting that the spectral function So{E) 
is equal to the particle (hole) spectral function at positive (negative) energies. 
The only differences from the results of Ref. [18] are the insertion of the average 
energy e and an energy shift in A^^{—u!, q) and in A^^J^—uj, q). 

The problem of defining a projection operator related to a single continuous 
eigenvalue e of the residual nucleus is of strictly mathematical nature. It is 
solved in the Appendix associating with e a set of approximate eigenvectors 
depending on an index 77 which describes the accuracy of the approximation 
(increasing for rj — > +0). They can replace the exact eigenvectors |e) in the 
whole formalism for two reasons, i) In the limit for 77 +0 they satisfy a 
completeness relation analogous to Eq. (80). ii) One can associate with e a 
projection operator since the approximate eigenvectors are normalizable. 

3.4-4 N on- energy-weighted sum rules 

We shall prove here that the approximation of Eq. (78), relating the spectral 
functions Sn{E) to (So(£^), has no effect on the sum rules, both for the inco- 
herent and the coherent contributions to the hadronic tensor. Since in Eq. (69) 
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we proved that the interference term W^{uj, q) does not contribute to the sum 
rule, Eqs. (81), (83), and (85) should satisfy the same sum rules as the exact 
hadronic tensor. 

Remembering Eq. (49) and using Eq. (13), i.e. 

/+0O 
dESo{E)^l, (89) 
-oo 

Eq. (81) yields 



/OO POO 
-0 J-oo 

/oo "I 
dcD5o(cD) =ReTr[Xq]. 
-oo J 



= ReTr 



By Eqs. (82) and (52), one has 



Ti\K^]= / dpdfc(^o|a;_qapaj,ak-q|^o) = (^o|J°^(g)J°(q)|^o) 



(90) 



(91) 



and so 

s°(q) = {MJ'\<DJ'm^^). (92) 

according to the exact sum rule of Eq. (48). Likewise, using Eqs. (83) and 
(84), one has 



/oo 
da; (a;, q) = Tr [K^] =Tr[K] = Z+ 1, 
-0 



(93) 



according to the exact sum rule of Eq. (49). It follows that also the sum rule 
for the coherent part is in accordance with Eq. (50). 



3.4-5 Practical approximations for the direct contribution to the hadronic 
tensor 

We have previously expressed the direct contribution W^{q) to the hadronic 
tensor as the sum of the incoherent and coherent parts 

Wr,,{u, q) = Al,{u, q) + Al,{-u, q), (94) 

W^^^iuj, q) = Al^{u, q) + A\{-uj, q). (95) 

In this Subsection, we profit from some results of Rcf. [18] to simplify W^{uj, q) 
in view of practical calculations. We notice that the terms of Eqs. (94) and 
(95) at positive (negative) arguments correspond to the particle (hole) contri- 
butions of the quoted reference. 
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In the calculations based on the Green's function approach, the coherent part 
of the hadronic tensor is usually disregarded since it is too difficult to evaluate. 
Consequently, only the sum rule EPj,(q) = Z + 1 of (93) is relevant. This is 
considered correct at high momentum transfers, schematically for q > kp. In 
Sec. 10 of Ref. [18], it was shown that this approximation is excessive even for 
uncorrelated systems, since only A^^^[u!,q) is negligible at g > A;^, whereas 
A^^Y^{—uj,q) is still sizable for < q < 2k-p. Fortunately, the latter term is 
largely cancelled by A^^^{—u}, q) and this produces the simplification 

W''{uj,q)^Al^{uj,q) at q>kp. (96) 

The energy shifts that are introduced here in the arguments of A?^^{—ui,q) 
and A^^-^{—uj, q) do not influence the cancellation found in Ref. [18]. 

Due to this cancellation, the sum rule for W^{uj, q) is deprived of the contri- 
bution of /lPp(— a;,q), which is positive as shown in Eq. (87). Thus, the sum 
rule for the total hadronic tensor, which coincides with that of W^{lo, q) since 
the interference term does not contribute, reduces to 

E(q) < Z + 1 for kp < q < Ikp. (97) 
3. 5 Inclusion of the interference term 

The interference term W^{u!,q), defined in Eqs. (65) and (67), gives no con- 
tribution in absence of final state correlations, since in this case S{E — Hn) 
does not connect the projection operators Pn and Q„. In the Green's function 
approach W^{uj, q) is usually disregarded also in presence of final-state corre- 
lations, since it does not seem reducible to a one-body expression. As already 
noticed, this approximation does not affect the sum rule. We give here an ap- 
proximated expression of W^iuj, q), which is reducible to a one-body expression 
and does not modify the sum rule. This result has a conceptual relevance and 
will be useful to introduce the empirical optical-model potential. 

Here, we are only interested in the interference term A^{uj,q) to be added 
to Eq. (96). Therefore, we start from the expression (67) of W^{u!,q) and 
disregard the contribution A^{—iJ,q). Then, we introduce an approximation 
which allows the reduction of A^{u,q) to a one-body expression and, finally, 
we retain only its incoherent part. 

We insert into the expression (65) of A^{u;, q) the Eq. (26), i.e., 

5{E - H,;) ^ ^[Gi{E) - G^{E)i (98) 
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with 

Gn{E) = ^ }. . (99) 
E- Hn + iT] 

Then, we use the approximated relations 

PnGn{E)QnJ%q)\iJo) = -PnGn{E)PnM'^{E)j\q)\^Po), (100) 

Pr,Gi{E)Q^j\q)\i;,) = -P^Gi{E)P^{M'J{E)j\q)\i;,), (101) 

where M^{E) is the energy derivative of the many-body self-energy of Eq. (33), 
obtaining 



1 f 
A\uj,q)^ H^e/ dp(V'olap-ql^) 



with a;„ = a; -|- £^0 — 

Equations (100) and (101) have the same structure as Eq. (59) of Ref. [28] and 
can be deduced with the same method replacing the Hamiltonian H by Hn 

and the Feshbach's projection operators by the extended ones. As remarked 
in the quoted reference, the Eqs. (100) and (101) must be used inside the 
matrix elements of Eq. (102) and hold in the region of the quasielastic peak 
at intermediate and high energies. 

Operating as in Sec. 3.4.1 and retaining only the incoherent part of Eq. (102), 
we obtain 



X /(q)|0n), (103) 



where M.n{E) is the one-body expression of the self-energy. Using the relation 

deduced from Eq. (11), and the relation [see Eq. (19)] 

Mn{E)^E-[g^{E)]-^-T, (105) 
Eq. (78) is extended to the relations 

^„(u; + ^0 - e„) ~ Qo{Cj), Mn{oj + Eq- tn) ^ Mo{u;) , 

u ^ u -\- Eq — eQ — e. (106) 
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Therefore, the contribution of the continuous spectrum can be included as in 
Sec. 3.4.3 to yield 



Al^{u;,q) = —ReTr 



TT 



2i 



Making exphcit the real part of Eq. (107), we have 



(107) 



'[gl{M'oV + {M'oVglm 



2i 

where M.'q{E) and QoiE) are symmetrically arranged. 
Equation (107) imphes 

'^^ dujAl^{uj,q) = 0, 



(108) 



(109) 



in keeping with Eq. (69). In fact, denoting by Cqo the large circle with center 
in the origin, one has 



+ 00 



dE[gl{E){M'oy{E)-go{E)M',{E)] = / dz go{z)M'o{z) = 0, (110) 

Coo 



where the first equality is due to the analyticity of go{z)M.'Q{z) in the complex 
plane, except for cuts and poles on the real axis, and the latter one follows 
from the fast decrease of go{z)A4o{z) at infinity (more than l/l-zj). 



3.6 Practical approximation for the total hadronic tensor 



In this Subsection we consider the total hadronic tensor as obtained by the sum 
of the direct contribution of Eq. (96) with the interference term of Eq. (108): 



W-(u;,g)=A°nc(^,9)+^L>,q)- 
Using in Eq. (87), the relation (see Eq. (12)) 



(111) 



(112) 



one has 



, Co ^ u; + Eq — €o — i, (113) 
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where 

g^siE) = [1 - M'.iE)]'/' g^m^ - M'oiE)]'/'. (114) 
We call ^eff(-^) "effective Greeen's function". 

Strictly, Eq. (114) should involve in a symmetrical way the factor 1—Mq{E)/2, 
instead of the square root operator. The replacement by the latter has been 
done to feature the more fundamental quantity 1 — A4'q{E), related to the 
effective a;-mass (sec Ref. [35]), and is justified by the fact that M.q{E) is 
small in the energy region of interest. For instance, in the scattering p— ^°Ca 
the approximation is very good at proton energies beyond 50 MeV (see Fig. 2 
of Ref. [27], where the result has been obtained using the dispersion relation 
specific of the self-energy). Concerning the precise meaning of the square root 
operator of Eq. (114) , we remark that it is trivial if M.q{E) is local (as supposed 
in many models) and that in any case the square root can be defined without 
problems as a power series for ||A1q(£^)|| < 1. 

The effective Green's function shows an interesting property. While the energy 
derivative of go{E) satisfies the relation 

g',{E) = -go{E)[i - M',{E)]go{E), (115) 

the derivative of ^eff (-£"), performed with the help of the previous equation 
and disregarding the contribution of AiQ{E), as allowed by the quasilinear 
behavior of Mo{E) (see again Fig. 2 of Ref. [27]), yields 

g^E) ^ [1 - M',{E)Y/' g'.mi - M',{E)Y/' = -glAE). im 

This is the typical relation satisfied by the Green's function of an energy inde- 
pendent Hamiltonian. This property really holds. In fact g^slE) is invertible, 
since [1 — A4q(£')]~^/^ can be defined as a power series and, therefore, one 
can define the related self-energy MesiE) and the Hamiltonian hes{E) as in 
Eq. (19): 

M,^{E) = KsiE) - r , /leff(^) = E- g-^{E). (117) 
Note that the energy derivative 

MU{E) = 1 + g:^{E)g',^{E)g:^{E) (iis) 

is approximately equal to zero due to Eq. (116). This means that the hadronic 
tensor can be reduced to a one-body expression by means of the Green's 
function associated to an effective self-energy, nearly energy independent. 

Equation (113) holds under the same conditions as Eqs. (96) and (108), i.e. 
at intermediate and high energies, near the quasielastic peak and for q> k-p. 
Now, we arrange Eq. (113) in order to feature the role of the scalar operator 
3^{Q'i f) — exp(iqf • r). This is included in Kq as 

{r\Kq\r') = exp{-iq ■ r)K{r,r') exp{iq ■ r') = {r\f\q)Kf{q)\r'). (119) 
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Permuting the factors under the trace symbol and using the property 

Tr pt] = W\0], 
Eq. (113) can equivalently be written as 



(120) 



W{u;,q) = Tr[i^j°t(q)5eff(c:.)j°(q)] = --lmTr[Kf\q)g,s{Co)f{q)] , 

TV 

CU = LU + Eq - 60 - 6, (121) 

where 

Sesioj) = ^^[^Isi^) - ^efr(a^)]- (122) 

So far, we have considered only the scalar component J°(q) of the current 
in the momentum representation. The extension to an arbitrary one-body 
operator 

0{q) = J drdr'alO{q;r,r')ar' (123) 

is performed on sight. Thus, Eq. (121) is generalized to all the diagonal com- 
ponents of the hadronic tensor, provided that only one-body currents are con- 
sidered: 



W>'>'{u,q)^Tt[Kf\q)SM^)fiQ)]---^^'^[Kr^iQ)Ges{^^^^^ 

TT 



u^u + Eo-eo-e, (124) 



3. 7 Energy shift prescriptions 



In the literature two types of energy shifts have been proposed to connect the 
spectral functions Sn{E). 

a) Kinetic energy prescription [14,15,17,18]: 

Sn{(^ + Eq- en) = So{lj + Eq - - Se^) , Se^ = en- Co, (125) 

which produces an energy shift equal to e^— cq in the expression of the hadronic 
tensor. This approximation keeps the value of the argument of the spectral 
function corresponding to the kinetic energy of the emitted proton. This ex- 
plains the name. 

b) Total energy prescription [18]: 

Sniou + Eo- e„) = Sq{u} + Eo-eo- 5e„) , 5e„ = 0, (126) 
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which yields no energy shift in the expression of the hadronic tensor. This 
approximation keeps the value of the total energy of the system given by the 
proton and the residual nucleus in the state |n). 

The same prescriptions are applied to Gn{E) and A4n{E) according to Eqs. (104), 
and (105). 

The total energy prescription has the merit of being simple, since it requires 
no shift. In contrast, the kinetic energy prescription leads to an expression 
of the hadronic tensor which involves a convolution of Sq{E) with the hole 
spectral function related to iV'o) [18]- 

The kinetic energy prescription is the one that has been adopted by most 
authors [Refs. [14,15,17]]. Moreover, it is closely related to an expression pro- 
posed for high energy transfers [Refs. [29-34]]. However, we decided for an 
intermediate choice, due to the reasons explained below. 

Using the relations (see Eqs. (12) and (18)) 

S.(E) = ±lgUE) - 6,(^)1 , S„(E) = (127) 

the spectral function Sn{E) is decomposed as 

S^{E) = S'iiE) + 5jr(i?) (128) 

with 



St{E) - ri gi{E + ir^)g^{E + zry), (129) 

S'^iE) = -Gi{E)M\{E)g^{E) , M\ = l^[Mn{E) - Mi{E)\. (130) 
The hmit for 77 —> is understood, as usual, in Eq. (129). 



The pioneering paper of Ref. [14] has shown that for £' > Sf^{E) yields the 
elastic contribution, due to the overlaps (7^|ar|^k,n)i where I'f/'icn) describes 
a proton of kinetic energy k"^ /2m = E scattered by the residual nucleus in 
the state \n). The asymptotic condition of incident plane wave imposes that 
Sf^{E) must be related to Sq(E) so as to preserve k and hence the argument 
E. This is in favour of Eq. (125) rather than Eq. (126). In contrast, S^^{E) 
gives the inelastic contribution, due to the overlaps {n\aT\ipk,m) , with m n, 
which do not contain the plane wave. Therefore, the requirement of preserving 
k is not so stringent in this case, and prevails the exigency of conserving in 
Eq. (130) the size of A4^(£^), which is mainly determined by the number of 
the inelastic channels which arc open when a proton of kinetic energy k'^ /2m is 
scattered by \n). This number depends on the total energy /c^/2m + e„. Hence, 
Eq. (125) does not correctly yield the thresholds of the inelastic processes, 
and Eq. (126) is favoured to relate M\{E) to M\{E). Since the real and the 
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imaginary parts of the self-energy fulfill a dispersion relation (see Sec. 5.4 of 
Ref. [36]), the total energy prescription must be applied to A4n{E) and hence 
to GniE) and S':{E). 

Due to the above reasons, we adopt the kinetic energy prescription for S^{E) 
and the total energy one for S^^{E), i.e., 

Snico + Eo- en) - Sl\uj + Eo-eo-e)+ S'^{cu + Eo - eo). (131) 

Here, e is taken equal to the average excitation energy of the residual nucleus 
given in terms of the energies e„ and the spectroscopic factors A„, as 

- ^ En Ue - 6o) ^^32) 

According to Eq. (131), the effective spectral function Sn,es{E), defined as in 

Eq. (122), and which differs from Sn{E) by corrections involving only A^'„(i^), 
is related to the ground state effective spectral function Ses{E) of Eq. (122) 
as 

SnM^ + Eo- en) = SfsiuJ + Eo-eo-e)+ S'^{uj + Eo - eo), (133) 
where 

Sf^{E)=vgUE + ir,)g,s{E + iv): (134) 

S':^{E) = -gUE)Ml^{E)g,s{E) , MUe) ^ i(>leff - -Miff), (135) 

and Aics{E) is the effective mass operator defined in Eq. (117). Therefore, 
Eq. (124) becomes 

W^'^(u;,q)^Tr{Kf^iq)[Sf^(u+Eo-eo-e)+S'^(u;+Eo-eo)]f(q)} (136) 



3.8 Hadronic tensor in terms of the local equivalent potential 



We have expressed the hadronic tensor in terms of the self-energy, rather than 
in terms of the Feshbach's potential as done in Ref. [17], since the latter has 
a quite complicated spatial nonlocality [36] and its kernel is not symmetric 
in the coordinate representation. These drawbacks make problematic to find 
a local equivalent potential to be compared with the empirical optical-model 
potential. In contrast, the self-energy is symmetric, has a simpler nonlocal 
structure, and is considered the theoretical mean field most closely related to 
the empirical potential [36]. 

Here, we want to show that the hadronic tensor can be expressed in terms 
of the Perey-Buck local potential, phase equivalent to the self-energy. To this 
extent, we investigate the connection between the Green's function go{E) of 
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the self-energy A4o{E) and that related to its local equivalent potential. We 
start with a simple model of A4q{E), that is widely used to calculate the 
self-energy by means of dispersion relations (see Sec. 4.4 of Ref. [35]). It takes 
into account the fact that the nonlocality in the self-energy is gathered mainly 
in its statical part, which is hermitian, whereas the non-hermitian dynamical 
part has a weaker nonlocal structure. Thus, one sets 

Mo{E) = U + V{E), (137) 

where T>{E) is non-hermitian and local, whereas the energy independent term 
U is hermitian and has the Frahn-Lemmer nonlocal structure 

U{r,r')^H{\r-r'\)u(^^\r + r'\^ . (138) 

The precise shape of the nonlocality form factor H is not relevant for the next 
equations. Setting 

F(q2) = J ds exp(-iq ■ s)H{s), (139) 
the Perey-Buck potential [37,38] is defined by the implicit equation 

Vl{E) = UF(^E-Vl{E))+V{E). (140) 

Every eigenvector \xl{E)) ofVh{E) is related to the corresponding eigenvector 
\x{E)) of Aio{E) by the phase conserving relation 

\XL{E))^f{E)\x{E)), (141) 

where the inverse Perey factor f{E) [39,38] multiplies by the function 

f{E; r) = [l + U{r)F'{E - V^{E- r))] (142) 

Higher-order corrections [38], involving surface terms, have been disregarded 
in Eqs. (140) and (142). Analogously, the Green's functions related to Vl and 
Mo 

Q^[E) = [E-T -Vi^{E)+ iri\-\ (143) 
g^{E) = [E-T- Mo{E) + ir]]-\ (144) 

are connected by 

g^{E) = f{E)G,{E)f{E). (145) 
The energy derivatives of the two sides of Eq. (140) yield 

1 - M'^{E) = f{E)[l - V^{E)]f{E). (146) 

This relation is analogous to the relation among the effective masses m*, m^, 
and m,^ in the infinite nuclear matter (Sec. 3.6 of [35]). 
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In Eqs. (114), (117), and (122), we have defined the effective Green's function 
of Qo and the linked Hamiltonian, self-energy, and spectral function as 



GME) = [1 - M'oiE)]'/' Gom-^ - M',{E)]'/' , 
h,s{E) = E - g;s\E) , 
MME)^Ks{E)-T, 

Sen{E)^^.[gUE)-g,^{E)]. (147) 

Likewise, the corresponding quantities related to g^iE) and Vl(-E') are defined 

as 

gLME) = [1 - v'Lm'^' Sl{e)[i - vmf' , 

h^ME) = E - g^l^iE) , 

VLME)=hLME)-T. 

SlME) = ^iiSlME) - gLME)]. (148) 
Inserting Eq. (146) in the first Eq. (147) and then using Eq. (145), one has 

ges{E)^g^ME), (149) 

and then 

KfiiE) = hi^ME) , Mes{E) = VlME) , SefiiE) = SlME)- (150) 

We emphasize that this invariance property, which will be very useful below, 
is the mere consequence of having introduced into the hadronic tensor the 
contribution of the interference between different channels. 

As a consequence of Eqs. (149) and (150), the expression (136) of the hadronic 
tensor reads 

W^^iuj, q) = Tr {Kj^^^iq) [stM^ + E,^e,~e) 

+ SlA^ + Eo-eo)\f{q)} (151) 

with 

StAE) = vgiME + iri)gLME + iri), (152) 
Sts(E) = -glME)VlME)Gi.ME) , (153) 
where Vl^csiE) is the antihermitian part of VL,cff(-E'). Using the first two 
Eqs. (148) and Eq. (143), one readily obtains h^es{E) in terms of hi^{E): 

h^ME) = [1 - VdE)]-^/'[h^{E) - V'^{E)E][1 - V^{E)]-'/'. (154) 
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Therefore, in principle, the hadronic tensor, as well as its elastic and inelastic 
parts, can be calculated in terms of quantities related to Vl(-E). 

3.9 Practical calculations 

The local potential Vl(-E'), phase equivalent to the self-energy, can be identified 
with the empirical optical-model potential. Therefore, the symbol Vl(-E') will 
denote in the following the empirical potential. 

The use of two different energy shifts to treat the elastic and inelastic parts 
of the hadronic tensor represents a complication. In order to overcome this 
difficulty, it is sufficient to find a simple way of calculating the elastic part. 

We only deal with positive values of E, since we consider only scattering states. 
Let |xL7cfr(-^)) denote the eigenvector of /iL gg(£'), related to the eigenvalue E, 
and satisfying the condition of incoming scattering wave. The right hand side 
of Eq. (152) can be handled to obtain 

Ke^iE) = E \xi:kE)){xi:kE)i (iss) 

where the sum refers to the understood degeneracy indices. The nontrivial 
proof of Eq. (155) can be found in Sec. 4.12 of Ref. [36]. We remark that this 
proof extends, without any modification, to a general Green's function and its 
self-energy. 

Equation (155) involves the eigenvectors [xi^csiE)) of h\^^g{E), which is a 
complicated Hamiltonian, but can be expressed in terms of the correspond- 
ing eigenvectors |xl~^(-^)) of hl{E). In fact, by a simple substitution in the 
eigenvalue equation for /i[gg(£^) and using Eq. (154), one checks that these 
eigenvectors are related by the phase-conserving relation: 

\xi:kE)) = [1 - m\E)f'\x[-\E)). (156) 

Therefore, Eq. (155) reads 

Ke^iE) = E[l - {yj{E)f'\xi-\E)){xi-\Em - m{E)]'/'. (157) 

From this equation, the elastic contribution to the hadronic tensor can be 
calculated with the same difficulty as for the calculation of an integrated 
single-proton knockout. In contrast, no similar way for directly calculating 
the inelastic contribution is available. In order to overcome this difficulty, the 
hadronic tensor can be written as 

W^'^iuj, q) = Wr{oo, q) + W^^{u, q) (158) 
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with 



^--lmTr:{Kf\q)[l - {V^){E)]'/'g^{E)[l - {V^)m'^'nQ)l (159) 

and 

q) = Tr {Kj'^HqWlME - e) - StAEWiQ)}: (160) 
where E = lu + Eq — eq. 

In the following, W^^{uj, q) will be calculated using Eq. (157), and W^^{uj, q) 
will be obtained from the method of Ref. [17], based on the spectral represen- 
tation of the s.p. Green's function. 



4 Spectral representation of the hadronic tensor 

In this Section we consider the spectral representation of the s.p. Green's 
function which allows practical calculations of the hadron tensor of Eq. (159). 
Due to the complex nature of Vl(-E'), the spectral representation of Q-l{E) 
involves a biorthogonal expansion in terms of the eigenfunctions of hi,{E) and 
h\^{E). We consider the incoming wave scattering solutions of the eigenvalue 
equations, i.e., 

(^-4(^))lxU(£^)) = 0, (161) 
{£-K{E))\x^}{E))^Q. (162) 

The choice of incoming wave solutions is not strictly necessary, but it is conve- 
nient in order to have a closer comparison with the treatment of the exclusive 
reactions, where the final states fulfill this asymptotic condition and are the 
eigenfunctions | XL_i;(-E)) of h\^{E). 

The eigenfunctions of Eqs. (161) and (162) satisfy the biorthogonality condi- 
tion 

{xtkE) \ 5<^d{E)) = 5 {£ - £') , (163) 
and, in absence of bound eigenstates, the completeness relation 

f dE\xi:kE)){xVe{E)hl. (164) 

Equations (163) and (164) are the mathematical basis for the biorthogonal 
expansions. The contribution of possible bound state solutions of Eqs. (161) 
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and (162) can be disregarded in Eq. (164) since their effect on tlic liadron 
tensor is negligible at the energy and momentum transfers considered in this 
paper. 

Using Eqs. (164) and (162), one obtains the spectral representation 

g^{E) = f dS I xi:kE))^^^{xVsiE) I . (165) 

Therefore, Wi^{Lj,q) (Eq. (159)) can be written, after expanding the ODM 
K in terms of the natural orbitals, as 



Wriuj,q) 



— Im 



dS 



u; + Eq — €o — £ + ir] 



(166) 



where 



Tr{£, E) = n.{u. I f\q)Jl-Vm \ xU(^)) 
x{xVsiE)\\ll-V^{E)f{q)\u.) . 



(167) 



The quantities \uy) and rin are the natural orbitals and the occupation num- 
bers, defined in Eq. (5), and come from the natural expansion of the ODM K 
of Eq. (159). The hmit for 77 — > +0, understood before the integral of Eq. (166), 
can be calculated exploiting the standard symbolic relation 



lim 



(168) 



r,-.+o E -£ + irj 

where V denotes the principal value of the integral. Therefore, Eq. (166) reads 



wriu;, g) = E 

TT Jo 



ReT,^^(cu + Eo-eo,uj + Eo- eo) 
1 



dS- 



UJ + Eq - Cq- £ 



ImTI^''{S,u; + Eo-eo) 



, (169) 



which separately involves the real and imaginary parts of T^'*. 

The contribution coming from Eq. (160) can be calculated using Eq. (157), 
i.e. as 
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Wr{u,q)^Y.n.[{u^f\q)^]l-Vm\xtl{E)) 



X {xVEm^Jl-ViiE)f{q)\u.)] , (170) 

with E = uj + Eo — eQ and E — E—e.ln this equation, only the optical potential 
wave functions of Eq. (161) appear everywhere, and the calculation is similar 
to that of the integral of the exclusive cross section, but for the different energy 
values which are involved and the presence of occupation numbers and natural 
orbitals. 

Some remarks on Eq. (169) are in order. Let us disregard the square root cor- 
rection, due to interference effects, and the minor contribution of the integral 
over the energy. One has 

V 

Let us compare Eq. (171) with the corresponding elastic contribution due to 
»Sl(-E), given by Eq. (157) without the square root corrections, i.e. 



In Eq. (172) the attenuation of the strength, mathematically due to the imag- 
inary part of Vl(-E), is related to the flux lost toward the inelastic channels. 
In the inclusive response this attenuation must be compensated by a corre- 
sponding gain due to the inelastic contribution to W^^iy)^ q) . In the description 
provided by Eq. (171), including both elastic and inelastic contributions, the 
attenuation of strength of the first factor in the square bracket is compen- 
sated by the second factor, where the imaginary part of Vi,{E) has the effect 
of increasing the strength. We want to stress here that in the Green's function 
approach it is just the imaginary part of the optical potential which accounts 
for the redistribution of the strength among different channels. 

The main difference of the present approach with the one of Ref. [17] is that 
now we are able to include in the initial states the effect of correlations. In 
Ref. [17], indeed, the initial states were calculated as the overlaps between the 
target nucleus and the residual nucleus described as a hole state, corresponding 
to the separation energy taken from phenomenology, or computed through an 
independent-particle model. Here, the initial states are described through a 
realistic ODM, including correlations. However, this goal is obtained using a 
constant separation energy, i.e. the one corresponding to the ground state of 
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the residual nucleus, and this produces an undue shift of the cross section. 
The contribution of Eq. (160) is added in order to compensate the shift. 



5 Realistic one-body density and correlations 

The main issue of this paper is, together with a new and completely anti- 
symmetrized presentation of the Green's function approach, to investigate the 
effect of nuclear correlations in the inclusive quasielastic electron scattering. 
Correlations are included in the ODM, which is expressed within the natural 
orbital representation [40] in the form 

K{r, r') = {r\K\r') = ^ n,ul{r')u,{r). (173) 

V 

The calculations presented in the next Section arc done and compared for 
different realistic density matrices which include the contribution of short- 
range and tensor correlations and are obtained within different approaches. 
In particular, we consider the Jastrow correlation method (JCM) [41], the 
correlated basis function (CBF) method [42] , and the Green's function method 
(GFM) [43]. 

In Ref. [41] the ODM is obtained within the JCM in its low-order approxima- 
tion (LOA). The JCM incorporates the nucleon-nucleon (NN) short-range cor- 
relations (SRC) starting from the ansatz for the wave-function of A fermions 
[44]: 

^\v,,v,,...,va) = {Ca)-''' n F(|r,-r, |)$^(ri,r2,...,r^), (174) 

l<i<3<A 

where Ca is a normalization constant and is a single Slater determinant 
wave function built from harmonic-oscillator (HO) s.p. wave functions which 
depend on the oscillator parameter aosc.-, having the same value for both pro- 
tons and neutrons. Only central correlations are included in the correlation 
factor F(r), which is state-independent and has a simple Gaussian form 

F(r) = l-exp(-/3V), (175) 

where the correlation parameter j3 determines the healing distance. The LOA 
[45] keeps all terms up to the second order in (F — 1) and the first order in 
(F^ — 1) in such a way that the normalization of the density matrices is en- 
sured order by order. Under the above assumptions analytical expressions for 
the ODM and corresponding natural orbitals have been obtained in [41]. The 
values of the parameters aosc. and /5 have been obtained [41] phenomenologi- 
cally by fitting the experimental elastic form factor data. Thus, in the present 
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calculations the following values of the parameters are used: Oosc. = 0.59 fm ^, 
(3 = 1.43 fm-i for ^^O and aosc. = 0.52 im-\ (3 = 1.21 fm"! for ^oCa, respec- 
tively. 

In the CBF method the trial A-particle wave function has the form 



A 

n F{xi,xj) 

i<j=l 



^^{x,,...,xa), (176) 



where Xi are particle coordinates which contain spatial, spin, and isospin vari- 
ables, i5 is a symmetrization operator, and is an uncorrelated (Slater de- 
terminant) wave function normalized to unity and describing a closed-shell 
spherical system. The correlation factor F is generally written as 

F{xi,xj) = ^hn{\ri-rj\)dn (177) 

n 

with basic two-nucleon operators On inducing central, spin-spin, tensor and 
spin-orbit correlations, either with or without isospin exchange. The ODM gen- 
erated by a CBF- type correlated wave function for ^^O has been constructed 
in [42]. The LOA, which keeps terms up to the first order of the function 
h{xi,Xj]Xi,x'j) = F{x'i,x'j)F{xi,Xj) — 1, has been used. The s.p. orbits enter- 
ing the Slater determinant are taken from a Hartree-Fock calculation with 
the Skyrme-III effective force. The correlation factor F[xi,X2) obtained in 
[46] by variational Monte Carlo calculations with Argonne 1^14 NN forces has 
been used. The two-nucleon correlation factors were restricted to the central, 
spin-isospin and tensor-isospin operators. 

For a nucleus like ^^O with J — ground state angular momentum the ODM 
can easily be separated into submatrices of a given orbital angular momentum 
I and total angular momentum j. Within the GFM [47,48] the ODM in mo- 
mentum representation can be evaluated from the imaginary part of the s.p. 
Green's function by integrating 

Kijih, k2) = r dE-lm{gij{h, E)), (178) 

J-00 TT 



where the energy variable E corresponds to the energy difference between the 
ground state of the A particle system and the energies of the states in the 
{A — l)-particle system (negative E with large absolute value correspond to 
high excitation energies of the residual system) and ep is the Fermi energy. The 
s.p. Green's function Qij is obtained from the solution of the Dyson equation 

gij{h,k2;E)^gi°\h,k2;E) + J dk^ J dk,gj^\ki,ks; E) 

X AJ^uiks, k^- E)gij{k4, k2; E), (179) 
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where Q^^^ refers to the Hartree-Fock propagator and AE^j represents contri- 
butions to the real and imaginary part of the irreducible self-energy, which 
go beyond the Hartree-Fock approximation of the nucleon self-energy used to 
derive The results for the ODM have been analyzed in [43] in terms of 
the natural orbitals and the occupation numbers in ^^O. Within the natural 
orbital representation they can be determined by diagonalizing the ODM of 
the correlated system. The numerical results from [43] show that the ODM 
can be described quite accurately in terms of four natural orbitals for each 
partial wave Ij. 



6 Results and discussion 

The model presented in this paper has been applied to evaluate the response 
functions of the inclusive quasielastic electron scattering for the nuclei ^^O and 
^°Ca. For the nucleus ^^O, no data are available for the separated longitudinal 

and transverse response functions. However, we can compare the results pro- 
duced by several realistic ODM, obtained within different correlation methods. 
For the nucleus "^°Ca, experimental data are available and we can be compare 
our results with them. 

Calculations are performed in a nonrelativistic approach including the con- 
tributions of both Wi (Eq. (169)) and W2 (Eq. (170)), but disregarding the 
second term in Wi. The evaluation of this term, which contains the princi- 
pal value, requires the integration over all the eigenfunctions of the contin- 
uum spectrum of the optical potential and represents a quite complicate task. 
Moreover, it would be useless in a nonrelativistic approach where the con- 
tribution of this term is very small [17]. In Wi and W2 the sum over all the 
natural orbitals Ui, is involved and partial occupation numbers riu are included 
in the model. All the results presented in the following are normalized to the 
number of nucleons in the nucleus. The s.p. energies e„ and the spectroscopic 
factors A„ we have used in the calculation of the W2 contribution are those 
corresponding to the overlap states obtained in Ref. [49] for JCM, in [42] for 
CBF, and in [50] for GFM, following the method of Ref. [51]. 

The other ingredients in the matrix elements of Eqs. (167) and (170) are the 
same as those taken in our previous application of the Green's function ap- 
proach to the inclusive electron scattering (Ref. [17]) and in our treatment 
of the exclusive one-nucleon knockout, which is based on the distorted wave 
impulse approximation and was widely and successfully applied to the anal- 
ysis of {e.e'p) data [1,52]. The one-body nuclear current operator is given 
by the nonrelativistic approximation of Ref. [53], x X ^-^e eigenfunctions 
of a phenomenological spin-dependent optical potential and of its hermitian 
conjugate, determined through a fit to elastic nucleon-nucleus scattering data 
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including cross sections and polarizations [54] . This allows a consistent treat- 
ment of FSI in the inclusive and in the exclusive electron scattering. The role 
of FSI in the Green's function approach was already investigated in Ref. [17] 
in a nonrelativistic framework and, more recently, in Ref. [11] in a relativistic 
framework and will not be discussed later on in this paper. 
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Fig. 1. Longitudinal (upper panel) and transverse (lower panel) response functions 
for the ^^0(e,e') reaction at q = 400 MeV/c. The curves indicated by JCM Wi 
and JCM Wi + W2 give the contributions of Wi and of the sum Wi + W2 (see 
Eqs. (169) and (170)) with the ODM of the JCM [41]. SM is obtained using the SM 
prescription of Ref. [17] with the phenomenological s.p. wave functions of Ref. [55]. 
SM Wi + W2 is the sum of the two terms Wi and W2 where the ODM is replaced 
by the SM prescription (Eq. (182)) with the phenomenological s.p. wave functions. 
SM HO Wi gives the contribution of Wi in the SM with harmonic oscillator wave 
functions. 
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A numerical example for the longitudinal and transverse response functions 
of the ^^0(e,e') reaction at g = 400 MeV/c is given in Fig. 1. The curves 
called JCM Wi and JCM Wi + W2 give the contribution of Wi and of the 
sum Wi + W2 and are calculated with the ODM within the JCM in its LOA, 
where only short-range correlations are included. The SM curve corresponds to 
the prescription adopted in our previous applications of the Green's function 
approach [17], where correlations are neglected and the components of the 
hadronic tensor are given by 

W^'^icu, q)^Yl ^^TZico + Eo-en,u; + Eo- e„), (180) 

n 

with 

TriE, E) = A„(v.„ I f\ q)sJl-Vm I xVEiE)) 

X {xi:kE) I - Vi.(^)/(g) I ^n) . (181) 

A pure shell model (SM) is assumed for the nuclear structure: the sum is over 
all the occupied states in the SM and a unitary spectral strength (A„ — 1) 
is taken for each s.p. state (/?„. In the calculations for the SM curve the wave 
functions are taken from a phenomenological Woods-Saxon potential [55] 
and e„ are the experimental excitation energies of the states n. 

The difference between the correlated JCM Wi + W2 and the uncorrelated SM 
results indicates that correlations give a redistribution of the strength (the 
total strength is conserved in all the calculations presented in this paper) and 
a reduction of the response functions by ~ 10% in the peak region. Part of 
the difference is however due to the different methods. If the sum Wi + W2 is 
calculated in a SM approach, where the ODM is replaced by 

{r\K\r') = J2K^*nir')Mr), (182) 

n 

the corresponding result SM Wi + W2, which is also displayed in the figure, 
lies between the SM and JCM Wi + W2 curves, and the reduction produced 
by SRC in the peak region is no more than ~ 5%. The different position of 
the peaks of the curves JCM Wi + W2 and SM is mainly due to the following 
reasons, i) In the JCM W1 + W2 curve the inelastic part of the spectral function 
Sn is calculated using the total energy prescription of Eq. (126) which does not 
produce a shift, differently from the kinetic energy prescription (125) adopted 
for the SM curve. The effect of these different prescriptions is indicated by 
the distance between the peaks of SM Wi + W2 and SM. ii) In the cases JCM 
Wi + W2 and SM the calculations use different values of e„. The resulting 
effect is indicated by the distance between the peaks of JCM Wi + W2 and 
SM Wi + W2. 
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The term W2 gives a shift of the calculated response functions that is im- 
portant to determine the position of the peak and to bring it close to the 
SM result. In the previous application of Ref. [17] the SM result was able to 
give a fair description of the size and shape of the experimental longitudinal 
response of the ^^C(e, e') reaction, in a range of momentum transfer between 
400 and 550 MeV/c. In contrast, the experimental transverse response was 
generally underestimated. For the transverse response, however, an important 
contribution might be given by two-body meson-exchange currents, which are 
not included in the present model based on the s.p. Green's function. 

The last curve drawn in Fig. 1 gives the contribution of Wi calculated in the 
SM approach (Eq. (182)) and with the HO single-particle wave functions used 
in the calculation of the ODM with the JCM (SM HO Wi). The comparison 
with the JCM Wi result indicates that correlation effects are still within ~ 
5%, but in this case they enhance the response functions in the peak region. 
Phenomenological Woods- Saxon wave functions certainly represent a more 
reliable basis for a SM calculation. A calculation with HO wave functions, 
however, may allow a more consistent comparison with the correlated ODM 
of the JCM and may give an idea of the uncertainty produced by different s.p. 
wave functions. 

The comparison between the uncorrelated SM and the correlated JCM results 
shows that the effects of correlations depend on the uncorrelated result that 
is considered for the comparison. The effects of SRC are however small and 
within the range of uncertainty produced by the choice of the theoretical ingre- 
dients. From this point of view our results are in substantial agreement with 
those of Ref. [10], where the effects of SRC on the response functions of the 
quasielastic electron scattering were investigated in a different model. In Fig. 1 
correlation effects have a similar behavior on the longitudinal and transverse 
responses, while in Ref. [10] correlations act differently on the two responses. 
This apparent discrepancy might be explained by the different models used 
here and in Ref. [10], by the small effects produced by correlations, and by 
the uncertainties related to the other ingredients of the calculations. 

The longitudinal response functions of the ^^0(e, e') reaction ai q — 400 
MeV/c obtained with the ODM within the CBF [42] and the GFM [43] are 
shown in Fig. 2. In both calculations the ODM includes both short-range and 
tensor correlations. The two density matrices, which are obtained within dif- 
ferent correlation methods, give similar results. In comparison with the SM 
curve, correlations produce a redistribution of the strength and a reduction 
of the response in the peak region by ~ 15%. In comparison with the SM 
Wi + W2 result (that is not shown Fig. 2, but is given in Fig. 1) the reduction 
is within 10%, but anyhow larger than for the JCM density matrix, where 
only SRC are included. The term W2 gives a shift of the response function 
that brings the position of the maximum toward the SM result. 
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Fig. 2. Longitudinal response functions for the ^^0(e, e') reaction at g = 400 MeV/c. 
The contributions of Wx and of the sum Wx + W2 with the ODM of the CBF [42] 
and of the GFM [43] are displayed in the upper and lower panels, respectively. The 
SM result is the same as in Fig. 1. 



The longitudinal response functions of the ^'^Ca(e, e') reaction at g = 450 
MeV/c is shown in Fig. 3. The JCM Wi and JCM Wi + results are com- 
pared with the SM and SM W\ + 14^2 ones and with the available data [56,2]. 
The main difference between the different results is in the position of the 
maximum. The contribution of shifts the response by ~ 15 MeV at larger 
values of the energy transfer. The difference between the SM W\ + W2 and SM 
results and between the JCM W\ + W2 and SM W\ + results indicates the 
sensitivity of the position of the peak to the different prescriptions used for 
the energy shifts and to the values of e„ used in the calculations. These effects 
are more important in a heavier nucleus like ^°Ca than in ^^O. In comparison 
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Fig. 3. Longitudinal response fTinctions for the '^^Ca{e, e') reaction at q = 450 
MeV/c. The contributions of Wi and of the sum Wi + W2 with the ODM of the 
JCM [41] are displayed together with the SM and SM Wi + W2 results produced 
by the phenomenological s.p. wave functions of Ref. [55]. The Saclay data (open 
squares) are from Ref. [56], the MIT-Bates (black circles) are from Ref. [2]. 

with data, our results are closer to the Saclay data. 



7 Summciry and conclusions 

In this paper we have presented a completely antisymmetrized Green's func- 
tion approach to inclusive quasielastic electron scattering. The main goals are 
the following. 

i) To express the final state interaction in terms of the self-energy, rather than 
in terms of the Feshbach optical potential as was done in previous papers, 
since the mass operator is more closely related to the empirical optical-model 
potentials. Therefore, wc have developed a theoretical approach based on ex- 
tended projection operators, as in Ref. [18]. 

ii) To include an approximated treatment of the interference between differ- 
ent reaction channels, which is usually disregarded. This has required some 
modifications of the treatment of Ref. [18], resulting in a different separation 
of the hadronic tensor into a direct and an interference part. The latter one 
is essential to explain the replacement of the self-energy with the empirical 
potential. 
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iii) To separate the elastic and inelastic contributions to the hadronic tensor, 
in order to introduce different approximations concerning the dependence on 
the state of the residual nucleus. 

iv) To include in the calculation the effects of correlations in the target nucleus 
through the use of a realistic one-body density matrix. 

The method has been applied to the reactions ^^0(e, e') and ^°Ca(e, e'). Results 
produced by density matrices including short-range as well as tensor correla- 
tions have been compared. The effects of short range correlations are small 
(within ~ 5%) and within the range of uncertainty related to the choice of 
other ingredients of the calculation (e.g. the s.p. wave functions). This result is 
in substantial agreement with previous investigations [10], but a similar effect 
is found in the present work on the longitudinal and the transverse response 
functions. Stronger effects are obtained when also tensor correlations are in- 
cluded. For the ^^0(e,e') reaction at g = 450 MeV/c the height of the peak 
is reduced by ~ 10% when both types of correlations are considered. A cor- 
rection dependent on the s.p. energies is necessary [the term W2 in Eq. (170)] 
in order to determine the position of the peak. The comparison of the present 
results with the available experimental data for the nucleus ^''Ca gives a better 
agreement, in the longitudinal response, with the Saclay data than with the 
MIT-Bates ones. 

In order to get deeper insight into the dependence of the inclusive electron 
scattering on correlations, one should also compute the contribution of two- 
body currents and their interference with tensor correlations, that seems to 
give a larger effect[12,9,13]. In order to accomplish this task, however, a new 
approach, based on the two-particle Green's function and including the two- 
body density matrix, should be developed. 



A APPENDIX 

A real number a belongs to the continuous spectrum of a self-adjoint operator 
A if and only if: 

i) it does not belong to the discrete spectrum, 

ii) for every 77 > 0, one can find a normalizable vector lipa,^) a-nd a constant 
C(j such that 

ii(«-^)K.)ir<c^aiiK.)nv (A.i) 

The vectors \ipa,r]) are called "approximate eigenvectors of A related to the 
eigenvalue a" [see Sec. 8.1 of Ref [57]]. Now, we construct a set of approximate 
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eigenvectors of the Hamiltonian H of the residual nucleus, which, added to 
the exact bound eigenvectors \n), satisfy a completeness relation in the limit 
f] +0. We remark that this property is not an automatic consequence of 
Eq. (A.l). 

Let Tic be the Hilbert subspace spanned by the non-normalizable eigenvec- 
tors |e) of the continuous spectrum of H which, for simplicity, is assumed 
to coincide with [0, +00). Let {|Mfc)} be a basis in Tic consisting of normal- 
izable vectors such that {e\uk) are functions of e bounded, continuous, and 
everywhere different from zero. This can be realized weighing the eigenvec- 
tors |e) by means of a complete orthonormal set of functions having the same 
properties. We set 



\e,k;7]) = J^e-H + t7])-^\uk). (A.2) 

V TT 

Note that the vectors \e,k;rj) and \n) are orthogonal. 

Theorem 1 . The vectors \e,k]r]) are approximate eigenvectors of H. 



Proof. One has 



lim II |e,fc; 77)11^=1 (eK)lVO (A.3) 



smce 



\\\e,k;ri)r = T de'L{e - e';vW\uk)\\ Hx;ri) ^ (A.4) 

and, in the distributional sense of the limit, 

hm^L{x;r))^5{x). (A.5) 



Moreover, the following inequality holds: 

1/ I (^-^ 

TT (e — Hy + r]^ tt tt 



e-H)\e,k;7j)f=%,\-^^-JP^\u,) < ^{u,\uk) = ^. (A.6) 



Therefore, one has 
and Eq. (A.l) is satisfied. □ 



H^Mi^^^^.O (A.7) 

ri-*+o |||e, A;;77)||2 ^ ^ 



Theorem 2 . In the limit for rj — > -|-0, one has the completeness relation 

POO 

+ lini V / de\e,k;ri){e,k;r]\ ^ 1, (A.8) 
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where the convergence is understood in the weak sense, i.e., inside a scalar 
product between normahzable states. 



Proof . Due to the orthogonahty between the vectors \e,k]r]) and it is 
sufficient to prove Eq. (A. 8) inside the scalar product between vectors |0) and 
\x) belonging to Tic, where both states {|e)} and {1^^)} are complete. Thus 
the sum over n does not contribute, and one considers only 

Y^j^ de{(j)\e,k;v){e,k;7]\x) = -J^ (0| _ ^ |x) 

poo roo 

= de de'L{e-e';r)Me'){e'\x) 

Jo Jo 

= de'{<P\e'){e'\x) dxL{x;7j), (A.9) 

Jo J-e' 

where in the last step we have exchanged the integrals, according to the Fu- 
bini theorem, since |L(e — e'; 77)(0|e')(e'|x)| is Lebesgue summable in e and e'. 
Observing that 

dxL{x-7])\< \ W){^\X) / dxL{x-v)\ 

-e' J —oo 

^W){e'\x). (A.IO) 

where the last term is Lebesgue summable and does not depend on 77, we can 
use the dominated convergence theorem to take the limit for t] +0 within 
the first integral of the last term in Eq. (A.9). Thus, using Eq. (A. 5), one 
obtains 



/•oo roo roo 

hm V / de(0|e,fc;r?)(e,A;;r?|x)= / de' (0|e')(e'|x) / dx8{x) 

r)^+0 ^ JQ Jo J-e' 

roo 

= / de'(0|e')(e'|x) = (0|x),V|0),|x)e7ic. □ (A.ll) 
Jo 



Substituting the completeness relation of Eq. (53) by (A. 8) and introducing 
the projection operators 

P{e,k;7^)= J dpap|6,fc;r;)^^(e,fc;r;|ap, |e,fc;r7)^ = (A.12) 

the contribution of the continuous spectrum is recovered adding to Eq. (64) 
the term 
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lim^^ J deRe J dpdp' (-^ohp-qk, ^; ^) 

^ k 

x^{e,k;r]\apS{uj + Eq- e- H^)ap>\e,k;r))^ {e,k;r)\ap>J^{q)\ipo), (A.13) 
where He is defined analogously to Hn (see Eq. (25)), i.e. as 

He^H-e in H^^^^^ , H^^e-H in n^^-^\ (A.14) 
Thus, using the approximation 

''{e,k-r]\a^5{uj + E^-e-H,)a^>\e,k-r))'' c^{p\S^{uj + Eo-eQ-e)\p'), (A.15) 
analogous to Eq. (78), and using again Eq. (A. 8) one recovers Eqs. (81)-(86). 
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